March, 2012 



(N 



£ 



OS 
(N 



> 



in 

en 



A note on perturbation series in 
supersymmetric gauge theories 

Jorge G. RussoZ 



Perimeter Institute for Theoretical Physics, Waterloo, Ontario, N2L 2Y5, Canada 
<N ■ 

Institucio Catalana de Recerca i Estudis Avangats (ICREA) 

Pg. Lluis Companys, 23, 08010 Barcelona, Spain 



Abstract 



Exact results in supersymmetric Chern-Simons and M = 2 Yang-Mills theories 
can be used to examine the quantum behavior of observables and the structure 
of the perturbative series. For the U(2) x U(2) ABJM model, we determine the 
asymptotic behavior of the perturbative series for the partition function and write 
it as a Borel transform. Similar results are obtained for M = 2 SU(2) super Yang- 
Mills theory with four fundamental flavors and in N = 2* super Yang-Mills theory, 



for the partition function as well as for the expectation values for Wilson loop and 't 
Hooft loop operators (in the and 1 instanton sectors). In all examples, one has an 
alternate perturbation series where the coefficient of the nth term increases as n\, 
and the perturbation series are Borel summable. We also calculate the expectation 
value for a Wilson loop operator in the M = 2* SU(N) theory at large N in different 
regimes of the 't Hooft gauge coupling and mass parameter. For large masses, the 
calculation reproduces the running gauge coupling for the pure J\f = 2 SYM theory. 

X 



*On leave of absence from Universitat de Barcelona and Institute of Cosmos Sciences, Barcelona, 
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1 Introduction 

An intriguing feature of quantum field theory is that the number of Feynman diagrams of 
nth order increases as n\ [HE]- F rom early times, this led to the belief that perturbation 
series is divergent. In some cases it may be rendered convergent by the extensively studied 
technique of Borel transformation [3j. If one has the series 



Z(g) = y £cn9 n , 

by inserting a "1" , 



ra=0 



1 = 1/ die'* t n , 



ni 



o 



and commuting integral and sum, one finds a new expression 

Z(g) = / dt e*B(tg) , B(x) = V -. x n , (1.1) 

where the expansion series for B(x) has better convergence properties. In particular, if 
c n ~ n\ then B(x) will have finite radius of convergence. In this case B(x) will have at 
least one singularity somewhere in the complex (Borel) plane at some finite x. As long as 
this singularity does not lie on the real positive axes, the integral over t can be performed, 
obtaining a finite expression for Z(g) out of the original divergent expression. The trick 
of course lies in the illegitimate step of commuting the infinite sum and the integration; 
although it may appear as an artificial trick, the examples considered in this paper give 
a new justification in support of the Borel resummation procedure. 

Some singularities of the Borel transform B(x) are associated with instanton solu- 
tions [2j. The special thing about instantons is that they tend to produce singularities 



on the real axes of B(x). In non-abelian gauge theories, instantons may produce singu- 
larities on the negative real axes, in such a case they are harmless. Configurations of 
instantons and anti-instantons of vanishing topological charge in current-current correla- 
tion functions are believed to produce singularities on the positive real axes. There can 
also be singularities which lie away from the real axes, and these typically correspond to 
(unphysical) complex field configurations. In non-abelian gauge theories the main prob- 
lem in applying the Borel transform arises from renormalons [3], associated with terms 
in the operator product expansion. They can produce singularities on the positive real 
axes, leading to an ill-defined Borel integral [J]. This obstruction has been extensively 
studied (see e.g. [5] and references therein) leading to the theory of resurgent functions 
and there are prescriptions to construct well-defined quantities by deforming the contour. 
This gives rise to an imaginary part which in some cases describes a physical instability 
of the system (an example is the ground state energy of the cubic anharmonic oscillator). 

In non-abelian gauge theories in four dimensions, many of these considerations are 
incomplete insofar as one does not know in general closed analytic expressions to all orders 
in the coupling for any process. In some cases, even obtaining partial results require a 
complicated analysis. However, the recent exact results in supersymmetric gauge theories 
obtained by localization techniques open the possibility of revisiting these old issues and 
understanding in a clean setting the origin of perturbation theory divergences. 

M = 2 supersymmetric gauge theories have a rich quantum dynamics which in many 
respects reproduces the expected physics of QCD. The exact low energy effective action 
constructed by Seiberg and Witten [6j[Tj shows chiral symmetry breaking and quark con- 
finement in the theory deformed by M = 1 preserving mass terms, via condensation of 
magnetic monopoles. More recently, exact results have been obtained for the partition 
function and for the expectation value of gauge invariant operators using localization [8,9J. 
Implementing localization requires the existence of a fermionic symmetry. One can show 
that the functional integral reduces to a submanifold of field configurations invariant un- 
der the fermionic symmetry that solve the saddle point equations. The complete effect of 
quantum fluctuations is then captured by the one-loop determinant. 

In M = 4 SYM the one-loop determinant is absent and as we will see this has a striking 
effect on the behavior of the perturbative series. One exact result in A^ = 4 SU(N) SYM 
is the expectation value of the circular Wilson loop operator obtained by Drukker and 
Gross [TU], extending the result of [H] beyond the planar limit. The result is expressed 
in terms of the following closed analytic expression: 

g 2 (l+iV) 

(W(C circlc )) = 6 " / dt e-H N ^h{Vtg) , (1.2) 

Nig Jo 

where I\ is the modified Bessel function. Comparing with (II. ip . we see that the integral 
is already given as a "Borel transform" , 

2e~ 



(W(C circlc )) = m J dt e-'Bitg 2 ) , (1.3) 



where 

B{x) = x N -^h{^) . (1.4) 

Consider the perturbation series of (W). Using the power series expansion of Bessel 
function J 1; we find 

<W(r )) e siv ( N + k)\g 2k 

J_f> 2n (-!)» " (-l)^(JV+fc)!(l + ^- fc 

Nl^ " 9 ' Cn 8™ ^ fc!(n-A;)!(Jfc+l)! 

n=0 fc=0 V J \ ' J 

For large n, c n ~ b n /n\, b = (N — 1)/8N. Thus the perturbative series for this observable 
has an infinite radius of convergence. As we will see this is not the case for gauge theories 
with less supersymmetry, where the perturbation series will be divergent. 

Consider now the planar N — y oo limit, with A = g 2 N fixed. Using (N + k)\/N\ m N k , 
we see that (W(C C i TC \ c )) —y 2/ 1 (vA)/vA, which shows that the series expansion in planar 
diagrams also has an infinite radius of convergence for this observable. 

In this paper we consider the following theories: the ABJM theory [12J, the M = 2 
superconformal SU(2) Yang-Mills theory with four fundamental hypermultiplets, and the 
M = 2* SU(N) theory obtained by adding a mass term for the adjoint hypermultiplet of 
M = 4 super Yang Mills (SYM) theory. Various properties of these M = 2 SYM theories 
are discussed in [7]. 

In section 2 we consider the exact expression for the partition function of U(2) x U(2) 
ABJM theory on S 3 obtained by Kapustin et al in [13]. We determine the asymptotic 
behavior of the perturbation series and show that it is summable in the sense of Borel. In 
section 3 we show that the partition function in the zero-instanton sector of the M = 2 
superconformal SU(2) Yang-Mills theory on S 4 can be written as a Borel transform. This 
formula explicitly exhibits the singularity structure and the convergence radius of the 
associated series. We also determine the asymptotic behavior of the perturbation series 
by a saddle point evaluation and discuss the one-instanton sector. The analysis is then 
extended to the expectation value of a circular Wilson loop operator. Section 4 contains 
several calculations on the M = 2* SYM theory. We first consider the SU(N) theory in the 
planar limit and compute the distribution of eigenvalues and the expectation value of the 
circular Wilson loop in different regimes of the gauge coupling and mass parameters. We 
then consider the SU{2) gauge theory. We discuss different properties for the partition 
function, Wilson loop and 't Hooft loop operators, we comment on the structure and 
asymptotic behavior of the perturbative series, the role of instantons and on the small 
and large mass behavior. Section 5 contains concluding remarks. 

2 ABJM matrix model 

1. The ABJM theory [12] is constructed by using two copies of Chern-Simons theory 
with gauge group U(N) x U(N) and opposite levels k, —k. The general formula for the 



partition function of ABJM theory on S 3 is given by 

ZABm - ^_ [ ^ a", IL^' (*?*)**'(?) eSE » , (,?-,?) (2 „ 

abjm 22NN]2 J (to)" (tor)" n,j ro sh 2 (^) ' l ' J 

1/k plays the role of coupling constant. By some simple manipulations the partition 
function can be casted in the following form [T^[T5] 

ZABJM = ^^ Yl i' 1 )^ / T^JAN- — " " ( Xl - Xa{l) \ ■ ( 2 ' 2 ) 

1 iy -*zS N J [Z7TK) riiCosh(f)cosh^^f^J 

Different aspects of the partition function have also been discussed in [TBIE] ■ 

We begin the discussion with a case where calculations can be made explicitly. The 
ABJM model with gauge group U{2) x U(2), i.e. N = 2. In this case we have 

Z ABJM (N = 2) = — J (2 ^ ) 2 cosh( ^ )cosh( ^ ) (^ " cogh 2 (xjgay 

' '1 - J(k)) , (2.3) 



32k 2 

where 



1 f 1 

J(k) = / dXidx 2 ; : ; : T7 r • (2-4) 

V ; 4W cosh m- cosh (a cosh 2 (*±=p) 



'C0sh(f ) COsh(f ) COsh 2 (2^) 

Using the Fourier integrals: 

2k 2 [ dt — J-— e it(xi ~ X2) 



cosh 2 (^g 51 ) J sinh(TrH) 

<&— r*^ = —^— (2.5) 

cosh ~ cosh(7rt) 



we can perform the integrals over Xi,x 2 . We get 



J(k) = -= / dv 9/ - . - . . , v = irtk . (2.6) 

tt 2 J cosh 2 (f ) sinh(t;) v ; 



This integral can be evaluated by residues, as done in [To]. The result is a finite sum of 
k different contributions of fc-depending terms. For the purpose of studying the power 
series expansion in 1/k, it is more convenient to work with the integral representation 
d2HD of J(k). We introduce 



7T 



2 i- dk .,. f, tanh(^ 



/(* )s -_/ F j(* ) = y*- jg; ^. (2.7) 



Using the formulas 

tanh (« ) = f ; 4-(4--l)^ F (^ 1 , (2.8) 

oo ? .2m-l i 

dv — — - = (-l) m+1 - vr 2m (4 m - l)5 2m , (2.9) 

.oo sinn (<v J m 

where £?2n are the Bernoulli numbers, we find 

^=-*|; i (- i )™^ (4m - 1)2 (T) 2m - (2 - io) 

Hence 

J(k) = -%I\k) = -% JT(-ir(2m - 1) (4- - I) 2 %i ( ^V m , (2.11) 
tt z tt z *-^ m[2m)\ V k I 

and 

vr 2 / 4tt 2 289vr 4 3844tt 6 \ 
64fc4 V 1 ~ 3P + 90fc 4 ~ 315fc 6 + '") ' ^ 2 ' 12 ^ 

It is interesting to examine the asymptotic behavior of this series. We use the familiar 
formula: 

lJ2n: - ^l) m+1 74((2m) 



(2m)! v ' (2tt] 



;-l) m+1 - ' (1 + 2~ 2m + 3- 2m + A~ 2m + ...) , (2.13) 

{2tt) 



We find the following behavior: 

(<-> \ 2m 
-H • (2.14) 

Thus the coefficient of l//c n increases as n\. The series (12.121) is therefore an asymptotic 
series in the sense of Poincare. As more terms of the series are included, the sum ap- 
proaches the true value of Zabjm(N = 2, k) until it begins to deviate leading to a divergent 
asymptotic behavior. Minimizing a generic term with respect to m one finds the optimal 
truncation (see discussion at the end of section 3). 

The series is divergent and at this point the main question is whether it is Borel 
summable. The associated Borel series can be defined in the standard way: one introduces 
a "1", 

' dt t 2m e-* , (2.15) 



(2m)! J 
6 



and write 



n=2 



(2.16) 
Because the series is divergent, commuting integral and sum does not give an equivalent 
result. Instead it defines the associated Borel series B(t), 

POO 

Z A Bm(N = 2,k)= dt e-* B(t/k) , (2.17) 

Jo 

m - (^ |(-l)"^|^ (4™ - l) 2 ^ i^r . (2.18, 

The perturbation series will be Borel summable if the integral ( 12. 17ft can be done, which 
requires that possible singularities in B(t) do not lie on the positive real t axes (or, if 
they do, they are integrable, like e.g. in the case of a Dirac delta function singularity). 
Asymptotically, B(t) behaves as follows 



m>l ^ ' 



2m 

(2.19) 



which has a finite radius of convergence. More precisely, inserting the full expansion (I2.13P , 
we find that B(t) is a sum of series with finite radii of convergence, the minimum radius 
being at R* — f so that all series are convergent for \t\ < R*. Other series have radius of 
convergence of the form R = ^, where p is a positive integer. There are singularities in 
Bit) on the imaginary axes at t = ±y and more generally at t = ^p where p G Z. 

Therefore we have shown that the perturbation series for the U(2) x U(2) ABJM 
partition function is Borel summable. This is an expected result, because there are no 
instantons in three dimensions. We recall that a singularity of B(t/k) on the real positive 
axes of t at t/k = a would have implied the presence of euclidean solutions with finite 
action e~* ~ e~ ak . 

2. We now consider an alternative treatment. In the present case, having closed analytic 
expressions, there is actually no need of artificially introducing a Borel associated sum by 
a Laplace transform. Indeed, one may define a slight modification of the Borel transform 
that works equally well and is more convenient in cases where the one-loop determinants 
involve hyperbolic trigonometric functions. Given an asymptotic series 

oo 

F(g) = J2f"9 n i ( 2 - 2 °) 

n=0 

one can introduce an associated series S(z) defined by 

f°° z 



with the expansion 

oo 

S(x) = J2*nX n , (2.22) 

n=0 

where s n is related to the original series coefficient f n by the formula 

/„ = 2~ n (2 n+2 - 1) C(n + 2)(n + 1)! s n . (2.23) 

Here we have used the formula 

/■OO ? .71+1 

2/ ^ — r ZT = 2- n (2" +2 -l)C(n + 2)(n + l)!. (2.24) 

J smh(f) 

Asymptotically, C(n) — >■ 1, therefore 

/„~(n+l)!s„, (2.25) 

like in the Borel case. Thus the transformation (12.211) has essentially the same effect 
of the standard Borel transform. In the present U{2) x £7(2) case, we do not know 
how to express Zabjm as a Borel transform in terms of closed analytic formulas, but we 
can express Z^bjm i n terms of the above 5-transform (12.211) . This has some obvious 
advantages. In particular, by (I2.2ip . one readily spots an associated series with finite 
radius of convergence, whose expression is determined by a closed analytic formula: 

z ABJM(JV = 2H _L_p„_£_ s( ^, (2 . 26) 

with 

S(x) = tanh 2 (x) . (2.27) 

By means of the sole structure of S one determines that: 

a) singularities of S arise at complex values for v , v = ik(2n + l)f with integer n. 

b) the series expansion of S in powers of 1/k has a finite radius of convergence ttt < f ■ 

c) the perturbation series for the partition function has an asymptotic behavior (2n)! for 
a term k~ 2n . 

d) the perturbation series is summable as the poles of the integrand do not occur at the 
positive real axes. 

3. Consider now the general U(N) x U(N) case. A way to understand the expected 
asymptotic behavior of the perturbative series is by a closer look at the partition function 
(12. ip . For this, it is convenient to write down the hyperbolic cosines in the familiar infinite 
product representation, 

00 / Ax 2 \ 

n=l x v ' / 



The perturbation series arises in expanding the integrand (except the exponential factor) 
in powers series of //j, v { and then integrating over /ij, z/j (again this involves an illegitimate 
commutation of sums and integrals). The series for the integrand have finite radii of 
convergence determined by the singularities at fa = Uj + i(2l + l)n where the cosh factors 
in the denominator vanish. Each integral over dfa or dvi of given factors /^ n \ v™ 1 times 
the exponential factor in ( 12. ip introduce factorials nj, mj which will lead to a divergent 
behavior of the perturbative series. The fact that the final divergent perturbative series is 
summable is obvious because it originates from the integral in ( 12. ip . which is well defined 
and convergent, in particular, it does not hit any of the singularities in the complex plane. 
A direct way to determine the asymptotic behavior of the perturbative series is by means 
of saddle-point techniques. Defining k = 1/k 2 , the coefficient d n of the n-th term n n of 
the power series in k is obtained by 

^ = ^-/^T^ABjM(iV, K ), (2.29) 

where the contour surrounds the point k = 0. For large n, the integrals over k and over 
fii, z/j are dominated by saddle points in the complex plane, e.g. /ij = Uj + i(2l + l)n + 
0(l/n), etc. We will employ this method in discussing the perturbative series for M = 2 
super Yang-Mills theories. 

3 J\f = 2 SU(2) theory with Nf = 4 massless hypermul- 
tiplets 

Consider first M = 2 SU(N) theory with A// = 2N fundamental massless flavors. This 
theory has superconformal symmetry. The partition function on S A is given by [8] 

Zscf = c^- (iY2 - 1} J ' d N - l a\{{a t - a f e'^^ Z 1 _ loop (a) \Z- mst {a, g 2 )\ 2 , (3.1) 



Kj 



where 



Z ^M= l< > u \' n , , H{x) = H(l + ?L) e"-, (3.2) 



YliH™^) ' v ' ii v ™ 



n=l 



and Y^i a i — 0. cn is a numerical constant independent of g. Following [21], we introduced 
the function H related to the Barnes G function by the relation 

H(x) = e- (1+7) * 2 G(l + ix)G(l - ix) , (3.3) 

where 7 is the Euler's constant. The instanton factor has a more complicated form and 
general formulas are given in [T9Tl2~T] . For the case of S 4 , one must fix the Nekrasov 
deformation parameters [19] to e\ = £2 = l/ r 5 where r is the radius of S 4 [22]. In what 
follows we set r = 1. 



Various properties of the partition function and circular Wilson loop for the large N 
SU(N) theory have been recently discussed in [231121] • In this section we focus on the 
SU(2) theory. For an SU{2) gauge group, the partition function reduces to 

4a 2 \2n 

9 " / »<> <>-< - I I 7: . ^L \ZZ (2) \ 2 , (3-4) 



/*00 2 ^^ / 1 

Z S ^ 2) = 1287T 5 / 2 <T 3 [ da aV^" ^ J] (1 ± gjL | 7- ^,2 

-°° n=l '- 1 - 



"SCF 



_i_ a z V 



where we have fixed the numerical constant cn to have -ZgcF^ = 1 at g = 0. Now 

tt (1 + ^) 2 " _ # 2 (2a 



n=l 



:i + s) 8n 



rSU(2) 

|G(l + i2a)| 4 
JPiay " |G(l + za)| 16 



(3.5) 



Consider the (k = 0)-instanton sector, obtained by setting |Z inst | — > 1 in (I3.4p . By 



5(7(2), 



the change of variable 
we get 

where 



a t 



a 



16vr 2 



7 SU(2) 
^SCF 



fc=0 



'7ra 



/•oo 

/ dt e~* B(at) 
Jo 



b(x) = v* n y ~ 



i + %f n 



n=l 



+ 



JL\Sn 



(3.6) 
(3.7) 

(3.8) 



Remarkably, in the zero instanton sector the partition function is already in the form of 
a Laplace transform of a Borel associated function. 

This compact form of B(x) explicitly exhibits the singularities in the Borel plane. 
There are singularities at 

at, = -n 2 , n = 1,2,... (3.9) 

The expansion in powers of at of B(at) has a finite radius of convergence given by |a£*| = 
1. Integrating over t adds a factor n\ to each n-th term a n of this expansion. Since 
the series for B(ott) has a finite radius of convergence, the n-th term of the perturbative 
series will have the asymptotic behavior n! (more precisely, (— l) n n\ since the singularity 
appears at a negative t). 

It is easy to compute the first few terms of the perturbative expansion. Using 



B(x) = yjx exp 



fc=2 



we find 

7 SU{2) 



x%(2k - 1) (2 2fe - 4) 



8505 



(3.10) 



k=0 



1 - 45a 2 C(3) + 525a 3 C(5) + a 4 (4((3) 2 - 7C(7)) 

8 

^a 5 (20C(3)C(5) - 17C(9)) 

135135 
16 



(3.11) 



-a 



(72C(3) 3 - 200C(5) 2 - 378C(3)C(7) + 341C(11)) + O (a 7 ) . 



10 



Note that no n factors appear in the perturbative expansion; it is an expansion in powers 
of a = g 2 /16n 2 and the associated Borel sum is convergent for \at\ < 1. The partition 
function is well defined for any a as there are no singularities in the integration region. 

The divergent perturbative series (13.111) is what one would have obtained if instead 
of using matrix model localization techniques the partition function had been computed 
using Feynman diagrams^ In the present context, it is clear that the divergence of the 
perturbative series emerges artificially, in the process of commuting the integral over t 
with the sum defining the expansion of B in powers of at. There is no divergence in the 
original integral defining the partition function. The singularities of B(at) lie on the neg- 
ative real axes, which means that the integral over t can be performed and therefore the 
divergent perturbative series computed by Feynman diagrams is summable in the Borel 
sense. It is the Feynman diagrammatic method which introduces the divergent behavior, 
because computing Feynman diagrams requires commuting the functional integral with 
the expansion in series of the exponential of the interaction terms. The convergent Borel 
associated sum is indeed the true result of the theory, as defined by the functional inte- 
gration; the divergent perturbative series here (and in the Feynman diagram calculation) 
arises artificially due to a mathematically illegitimate manipulation. 

Let us now compute the asymptotic behavior of the perturbative series by a saddle 
point calculation at large n. Writing 



we have 

1 f da 
a, 



'■n 



ySU(2) 



7 SU{2) 



= Y, d nOc n , (3.12) 

k=0 n=0 



oo 

5 



^- j ^jt-sym ■ - — -= r da I dte , ( 3 - 13 ) 

S(a,t) = t + (n + l)\oga - hgB(at) . (3.14) 

For large n, the saddle point equations in a and t are 

n tB'(at) aB'(at) 



a B(at) ' "~ B(at) 



(3.15) 



i.e. t* = n and a satisfying the equation aB' (an) / B(an) = 1. This equation has an 
infinite number of solutions, with an near the zeroes and singularities of B(an). The 
dominant saddle point (with smallest 5* and large eigenvalues of the Jacobian matrix), 
appears near the first pole, a*n = — 1 + 0{l/n). Computing 5* at this point, we find 

d n ~ e- 5 ^*'**) ~ (-l) n n n e~ n ~ (-1)™ n\ , (3.16) 

as anticipated above. 



^or example, the three-loop term with £(3) coefficient in the similar perturbative series for the Wilson 
loop, see (3.26), has been reproduced by computing three loop Feynman diagrams in |25) . 
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It is well known that the existence of instantons can lead to singularities in the Borel 
transform. Having closed analytic expressions, it is therefore interesting to study how the 
perturbation series behave in the non-zero instanton sectors and to determine the extent 
to which possible singularities produced by instantons could affect the convergence of inte- 
grals. Consider for concreteness the one-instanton sector. The one-instanton contribution 
is given by 

7 su(2) , -*£■ f dy ( y* \ 

J 2vr V [(y - a) 2 + 1) {(y + a) 2 + 1) / 



8^ 



e ~F~ (a 2 -3) . (3.18) 



Thus 



7 SU(2) 



oo 



.:> --" ' ' , :• .» ,2 H (2°) 



daa 2 e ~ (a 2 -3) 2 rra \ / . (3.19) 

k =i j-oo H*(a) 

This integral is well defined and convergent. Now the partition function is no longer in the 

form of a Laplace transform of a Borel function, because of the presence of the instanton 

weight e~~ . To write it as a Borel transform one would need to put it in the form 

7 SU{2) 



dz e" F(z) . (3.20) 

fe=l Jo 



This can be done by a change of coordinate, a 2 + 1 = z. We find 



F(z) = 6{z - 1) Vz^T (z - 4) 2 f(z - 1) . (3.21) 

This has a branch cut singularity at z = 1 (more generally, z = \k\ in the \k\ instanton 
sector), but obviously it does not affect the convergence of the integral. Let us now 
see how this singularity is manifested in the saddle-point evaluation. The analysis is as 
follows. In the 1-instanton sector, instead of (13. 14ft we would have 

S(a, a) = — + n log a - 2 log T }.. a \ + log a 2 - \og(a 2 - 3) 2 . (3.22) 

a H 4 (o) a 

The saddle points which lie near the singularities of log H (2a) /H 4 (a) are not significantly 
modified. The main difference is that now a new saddle point emerges, located at a 2 = 
a* = 1/n. Since a 2 ~ 0, terms log H (2a) /H A (a) and log(a 2 — 3) 2 are irrelevant. This is 
indeed the dominant saddle point. The corresponding saddle point in the \k\ instanton 
sector is located at a 2 = a* = \k\ln. One finds 



n" 



d n ~ — e- n ~n\ \k\~ n , (3.23) 

\k\ n 

which shows that the Borel function must have a singularity at z — \k\. This is precisely 
the location of the branch-cut singularity. 

In conclusion, the integral defining the partition function is convergent in the one- 
instanton sector. The instanton produces the expected singularity in the Borel function, 
which in the present case turns out to be harmless. 

12 



Wilson loop 

Similar considerations can be made for physical observables such as expectation values of 
Wilson loop operators. Do they exhibit the same asymptotic behavior? We restrict to 
the zero-instanton sector. For a circular Wilson loop located at the equator of S 4 in the 
spin-j representation, one has [8] 



oo ,„ 2 „ °° t\ i ia z \2n 



, f s i f°° 16 ^ 2 o 2 TT 1 + ^T) / V^ 

(tr i Pexp( / Adx + i%ds)) = Z' 1 / da a 2 e ~^~ \\ S^T Yl 



Air la 



n=l \ n z 



Consider the term e 47rfa (£ can be integer or half- integer) . We have 



(3.24) 



it.',.-, Jo°° dt e ' B(at) cosh(47r£va^) 



< e > = ^^ roc ,; .; W( \ ■ (3.25) 

J at e t B(at) 

This can be expanded in powers of a as done in [8 J . We find 

(e 4rfa ) = i + i27rW 2 + 40vr 4 a 2 £ 4 + « 3 (^— -7207r 2 / 2 C(3) 

+ a 4 (96vr 8 £ 8 - 5760vr 4 £ 4 C(3) + 12600vr 2 £ 2 C(5)) + 0(a 5 ) (3.26) 



This expression has been checked up to the three-loop a 3 order in [23] by explicit evalu- 
ation of the Feynman diagrams in perturbation theory. 

A term with a given power p of a is multiplied by a polynomial in £ of degree 2p. 
Terms with lower power of £ have coefficients which increase more rapidly with p. This 
is because each power £ 2s t s is accompanied by a factor l/(2s)! coming from the power 
expansion of cosh(47r£v / at). Thus the terms in the perturbative series for the Wilson loop 
that give rise to the most divergent behavior of the series are those with lower power of 
£ in the polynomial of degree 2p multiplying a p . For them, s — 1, so such terms clearly 
present the same asymptotic behavior (— l) n n\ that appeared in Z SC p . 

A different behavior can be found in a limit where £ — > oo and a — > with yfai fixed. 
This gives 



Airla\ . JO 



(e 4 " £a ) 



J °° dt e * \/aicosh(4:7i£\/ai) 



to 

4ir 2 a 



J °° dt e~* \fod 



(1 + 87r 2 a£ 2 ) 



a 2 £ 2 ^ 
1 + 2— I . (3.27) 



In this limit the one-loop determinant was set to 1, B(at) — > \fod. The result must then 
be the same as in TV = 4 SYM as described in the introduction. Indeed, by explicitly 
computing the integral (jl.2p with iV = 2, we exactly reproduce (I3.27P with £ = 1/2. 
Now the expansion in powers of a has an infinite radius of convergence and (e 47rfa ) is an 
exponentially increasing function of the parameter a£ 2 which starts at 1. 
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An interesting question regards the best resolution that the asymptotic perturbative 
series defining the Wilson loop can provide. This is done in the standard way by consid- 
ering the partial sum 



"0 



W no = J2 c na n , c n ~ {-l) n n\R- n . (3.28) 



n=l 



The optimal truncation is then obtained by using the Stirling approximation and minimiz- 
ing \c n \ ~ n\R~ n with respect to n. One finds no = R/a, for which the best resolution (or 
"non-perturbative ambiguity") is e = ^2nR/a e _iJ//a , where we used the Stirling formula. 
In the present case, we have closed expressions which allow an accurate determination of 
R, which is nothing but the radius of convergence of the associated Borel series. This 
gives R — 1, so that uq = IQ-k 2 / g 2 and the best resolution is e = y^2n/a e~ l l a ~ e _167r I 9 . 
For example, if g — 1, then a ~ 0.006; this gives n = 158 and e = 10 -67 . But if g = 8, 
a ~ 0.4, then n ~ 2.4 and the asymptotic formulas are not applicable. It should be 
noted that here the "non-perturbative ambiguity" e = ^2nR/a e~ R l a is unrelated to 
instanton contributions that would provide a "completion" of the series, since the origin 
of the divergence in the asymptotic series are singularities in the one- loop determinant. In 
particular, in the ABJM model discussed in section 2 there is a similar non-perturbative 
ambiguity e = 0(e~~) despite the fact that there are no instantons. 



4 M = 2* SU{N) Super Yang-Mills theory 

It is of interest to consider the behavior of the perturbative series and of loop operators in 
a case where there is no superconformal invariance. Typically, in such cases the one-loop 
determinant contains a divergence which needs to be renormalized [8]. One case where 
this divergence is innocuous is in the M = 2* SU(N) SYM, obtained by deformation of 
M = 4 by adding a mass M for the adjoint hypermultiplet. In this case the divergence is 
a numerical factor (in particular, it has no incidence on the expectation value of physical 
observables) . We remove this divergent numerical factor from the partition function by 

2 At ^ 

inserting e~~^~ in the infinite product over n in the integrand. The partition function on 
S 4 reads [8] 



Jd N ~'a Y[(a t - arfe-^^Z^^M) |Z^ (iV) | 2 , (4.1) 



?N=2 

= I u a 

i<j 



with 

H 2 (ai — dj) 



(1 - "'~" ; 



i<j 3 



nn 



(ai-a.,) 2 ^2n 



i< jn =i\\ - h (^- M H ( 1 + ^+Mh ) e ~ 
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M = 2* SU(N) super Yang-Mills theory interpolates between N = 4 SU(N) SYM and 
pure Af = 2 SU(N) SYM as the mass parameter is varied from zero to infinity. The 
interpolation seems to be smooth (for example, on R 4 the Seiberg-Witten curve presents 
no discontinuous or singular behavior for any M [7]). For large M, the M = 2* theory 
may be viewed as a regulated version of pure M = 2 theory where one has an UV cutoff 
determined by the mass. More generally, given an asymptotically free M = 2 theory, 
by adding a sufficient number of massive hyper mult iplets one gets a theory that is well 
defined in the UV and for large masses it provides a regulated version of the M = 2 
theory that one wishes to study. This interpretation was exploited in [8] to recover the j3 
function of M = 2 SYM from the partition function. This observation will be reproduced 
here in the planar limit. In addition, we will also study the small mass behavior and loop 
operators in the SU(2) theory. 

4.1 Large N limit and Wilson loop 

We consider planar limit, i.e. N — > oo with A = g 2 N fixed. In this limit \Z inst \ — > 1. 
We will compute the integral in the saddle point approximation. We use the notation 
of [21], where the similar calculation was carried out for the M = 2 superconformal 
theory discussed in the previous section. Our approach to solve the system is different: 
now there is an extra parameter, the mass, in addition to the 't Hooft coupling A = g 2 N, 
which will lead us to consider different regimes. 
The effective action for the eigenvalues is 



o 1 \-^ / , ,o , H (a,- 



a 



3/ 



H(cLi - CLj - M) H(di - CLj + M) 

(4.3) 



A N 



where 



i i<j 

The saddle-point equations are therefore 

S " "' h Yl \—z — K ^ ai ~ a i) + \ i K ^ ~ a i - M ) + K ( a i - % + M)) ) = , 

V \ Ulrj CI) A £l I 

(4.4) 

H'(x) 
K ( x ) = — FT^T =x{i/)(l + ix) +ip{l -ix) +2 7 ) 
n[x) 

oo 

= -2^(-l) fc C(2fc + l)x 2fc+1 , (4.5) 

fc=i 

and ip is the logarithmic derivative of the V function. 
Introducing as usual the eigenvalue density 

p{x) = —^5{x - ai) , (4.6) 
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the saddle-point equation can be written as the following integral equation: 

M / 1 1 \ 8tt 2 

dyp(y) K{x -y) + ~(K(x -y + M)+K(x-y-M)))= — x. (4.7) 

-n \x-y 2 J X 

Following [23], we do the transform 

^ dx 1 



/< 



z-x yjtf - X 2 



X 



The integral equation then takes the form 



p{x) = V, 



fl 2 — X 2 I — - 



In 1 f dy 



X 2n 2 J x-y^tf-y* 
x / dzp(z) (2K(y -z)- K{y - z + M) - K{y - z - M)) J . (4.1 



We will treat the second term on the RHS of ( 14. 8 p as a perturbation, which allows us 
to replace p(z) by no\/p 2 — z 2 where no = 2/(it[j, 2 ) is the normalization factor. The 
equation can be explicitly solved in two different regimes, p ^C M or /i ^> M. Since in 
the perturbative approach jjl ~ v^A, the first case will lead to a small A expansion, and 
the second case to a strong A expansion. We will now consider these two cases separately. 
By rescaling x,y,z — > px,py,pz, one gets new tilde variables which are of order 1. 
For |i « M we can expand K(u,(y — z) ± M) in powers of (y — z). One has 

2K(y - z) - K(y - z + M) - K(y - z - M) = -2K'(M) (y -z)+0((y- zf) . (4.9) 

By computing the integrals we find 

p(x) w y/p 2 - x 2 (y - - K'{M)\ . (4.10) 

The perturbative treatment of the saddle point equation is therefore justified for 

XK'(M) 



8tt 2 
The normalization condition gives 



< 1 . (4.11) 



A«« 1 = (4-12) 

27T,/l-^i 



The eigenvalue density (I4.10p can be written in the form 



p ( x )~^L v ^3^ ) (4.13) 
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with 

£ - t " i A "< M > ■ < 414 > 

K'(M) is a monotonically increasing function with the behavior 

K>(M) « ( 6C(3 > M2 r 10 «f /4 • M<<1 ' (4.15) 

\21ogMe 1 ^- 5 i I , M»l, 

where we have used (14.51) and the asymptotic form of the Barnes G function. 

G(l + iM)G(l - zM)~-^expf-(M 2 + i) log M+^- + 0(M~ 2 )) . (4.16) 

A is the Glaisher-Kinkelin constant, A w 1.28. Therefore, the condition /x <C M together 
with (14.111) require A ^C 1 for any M. For large M, one has a stronger requirement 
A<4vr 2 /logM. 

Defining g 2 R = \ c s/N and g 2 are = A /AT, for large M, one has 

47T 2 47T 2 AT 47T 2 

- —K'(Mr) « -^— - Anog(Mr) . (4.17) 



#R #bare ^ #ba 



where we have restored the dependence on the radius r of S 14 . Thus, in the setup M = 2* 
SYM is interpreted as M = 2 SYM without matter multiplets with a UV regulator M, g^ 
represents the renormalized gauge coupling. This interpretation leads to the familiar ex- 
pression for the one-loop /3 function. By replacing Mr by Vq/v, g]i{v) in (I4.17P represents 
the running gauge coupling of pure Af = 2 SYM, defined at an energy scale v. 

Returning to the M = 2* SYM theory for its own sake, we see from (14. 12j) that as M is 
gradually increased at fixed A<1, the interval (—//,//) of the eigenvalue distribution will 
expand. The approximation breaks down long before the denominator becomes singular. 

We can now use the above p(x) to calculate the expectation value of the Wilson loop 
operator at weak coupling. We find 



W(C cirdc ) = dx p(x) e 

J —n 



2nx 

mx pyx) e 
a* 

2/ 1 ( v /Ao»(A0) 



1 + » + »=ffl + i| + 0(V), (4.18) 



or simply 



±+\ 2 ( 

8 V 64tt 2 ' 192 

W(C M .) * 1 H- ^fl 4- &*£> 4- O (Ai) . (4.19) 

Thus in this approximation the complete M dependence in the Wilson loop is in X e g(M). 
It would be interesting to reproduce the non-trivial A 2 term proportional to K'(M) by a 
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perturbative gauge theory calculation. Note that the series ( 14 . 1 9 1) has an infinite radius 
of convergence, since it arises from the power series expansion of the Bessel function. The 
approximation of course breaks down for A ~ Ao = 8ir 2 /K'(M), so convergence is only 
demonstrated at A <C A . 

Using the above expansions of K'(M), one can explicitly see the small and large M 
behavior of the Wilson loop: 



W(C circlc ) 
W(C ciTclc ) 



M<<1 ~ 1+ 8 +A l^2^ + 192j +0 ^) ' 
^ 1+ A 2 /MM) 7 + 1 M 
8 V 327r 32vr 192 / 



A/>1 

where the second equation applies as long as AlogM <C 1 due to the condition (14. lip . 

Let us now consider the regime \i ^> M. In this case we can write 

2K(y -z)- K(y - z + M) - K(y - z - M) = - M 2 K"(y - z) + O (M 3 ) . (4.21) 

Using the rescaled variables y, z, it becomes clear that if A 3> 1 the integrals will get most 
of the contributions from the regions where the argument of K" is large. Therefore we 
can approximate K(x) ~ 2xlogx so that K"(x) ~ 2/x. Computing the integrals we now 
find 

p{x) « vV - x 2 T" (! - M 2 ) ■ ( 4 - 22 ) 

A 

Therefore the perturbative approach to the saddle point equation is justified for M 2 <C 1. 

Now we get 

fi « -^ . (4.23) 

27TV1 - M 2 V ; 

\i 3> M requires A ^> M 2 . Thus we have a calculation of the circular Wilson loop in the 
M = 2* theory for A ^> 1 at small M. In this regime, the circular Wilson loop in the 
M = 2* theory exhibits the following behavior 



^(C circle ) « e 2 ^ « exp 



2 



^(1 + — ) 



(4.24) 



This shows that the string tension T increases as a small mass parameter is turned on in 
the M = 4 theory, we get T « ^J(l + ^f ) at M < 1. 

4.2 A/" = 2* 577(2) theory 

A number of interesting properties can be derived explicitly when the gauge group is 
SU(2), For N = 2 the partition function reads 

+ n 2 J \ L ■+" n2 J e ™ 

(4.25) 
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Like in the M = 2 superconformal case, we note that in the zero-instanton sector the 
partition function can be written as a Laplace transform of the Borel associated function: 

we get 



7 Af=2* 



fc=0 



dt e _t B 



jAf=2* , 



at) , 



(4.26) 



where 



B 



Af=2* , 



X) = 



H 2 (2^) 



H(2y/x~ - M)H(2y/x~ + M) 

(1 + %) 2 ™ 



n 



i((i -*f) 2 + »)"((i+ >f) 2 + 



lx 



™ 2M 2 

e n 



and again we have defined 



ct £ , 



a = 



16tt 2 



(4.27) 



(4.28) 



The associated Borel series defined by the expansion of B N 2 (x) in powers of x has 



a finite radius of convergence determined by the first singularities that arise in B 
ai x = at = — t(1 ± iM) 2 . The associated Borel series is therefore convergent if 



■Af=2* 



X) 



N < |(i 



M 2 



Also in this case the perturbative series has asymptotic behavior n\, as follows from the 
fact that the series expansion of B(x) around x = has a finite radius of convergence. 
Because the singularities lie on the complex t plane, they do not affect the integration 
over the real positive axes of t. Therefore the integration is well defined (the integral is 
convergent) and the divergent perturbative series that one would obtain by computing 
Feynman diagrams is Borel summable. We have reproduced the n\ asymptotic behavior 
by a saddle-point calculation, similarly as done above in the M = 2 superconformal case. 
A systematic way to calculate the perturbative series is by writing B^ =T [x] as follows 



where 



B N=T { 

Pi(M) 

Pk(M) 



x) 



1 



H 2 (M) 



exp 



E 

fc=i 



-iy 



(Ax) k Pk (M) 



(4.29) 



K'{M) 
2((2k- 



D- 



1) 



tf(2fc-l)( M ) ^ k>X 



(4.30) 



(2ife-l)! 

This form is also suitable in order to take the large M limit, that we shall consider in 
discussing the Wilson and 't Hooft loop expectation values. For the first few terms one 
finds 

1 



7 N=2* 



H 2 (M) 



(1 + 6pi a + 30(p 2 - p 2 ) a 2 + 140(p? - 3 Pl p 2 + 2p 3 ) a 3 + 0(a 4 )) . 

(4.31) 
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Consider now the effects of instantons. In particular, the one-instanton contribution 
is given by [8] 

_s^£ 2(M - l) 2 (M 2 + a 2 - 2M + A) 

Zi-tort = e IP- -i ^- 1 , M = tM + l } (4.32) 

(4 + a 2 ) 

where we have set the radius of S 4 to 1. In the case M = one has Zi_i nst = and one 
recovers the result of M = 4 theory. 

Just as in the superconformal SYM case, in the one-instanton sector a new saddle point 
appears in the calculation of the coefficient d n of the rath term d n a n of the perturbative 
series, but it does not affect the convergence of the integral. We also see that the prefactor 
multiplying the instanton weight does not introduce singularities on the integration region 
which would otherwise represent an obstruction to computing the matrix integral. Indeed, 
the contribution (|4.32[) has a singularity in the complex plane at a = ±2i (and on the 
negative real axes for the integration variable t = a 2 /a). Therefore it is innocuous. 

4.2.1 Wilson loop 

Let us now consider the Wilson loop in the zero instanton sector: 

/OO 2 

da a 2 e'^MB^ 2 * (a 2 , M) cosh(2vrga) , (4.33) 

■OO 

where we have defined (cf. (14.29)) ) 

£ [ -^(2af p k [M) 



^" (ai - M) = wW)^ 



k=2 



(4.34) 



1 - 4pi(M) = - - AK'(M) . (4.35) 



a e g(M) a ' a 

This definition of a e g of course coincides with (14.141) when N = 2 (it extends to any N) . 
In perturbation theory we obtain 

(e 2 ^) = 1 + 37i 2 a cS q 2 + a 2 J-^f- + ^c^V (7vr 4 g 4 - 720p 2 ) + O (« c 4 ff ) . (4.36) 

z o 

which exhibits the dependence on the mass parameter M through a e a(M) and p2(M). 

One can also consider an expansion where a is small but q 2 (M) = q 2 a e e(M) is an 
arbitrary fixed number. We get 

(e 2 * qa ) = e^ (1 + 27r 2 g 2 - 8n 2 a 2 ca (M)p 2 (M)q 2 (2vr 4 g 4 + 15vr 2 g 2 + 15) + O (o^)) . 

(4.37) 
The leading term is once again the expectation value of the circular Wilson loop operator 
in M = 4 SYM since in the strict a — > limit with fixed q* the one-loop determinant is 
set to 1. 
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Let us apply these expressions to the study of the small and large M behavior. At 
small M, one has p k (M) « c k M 2 where c k = 2k{2k + l)({2k + 1). We find 

(e 2 ^ a ) = e^* (1 + 2it 2 ql - 1607r 2 C(5)a 2 ff M 2 g 2 (27r 4 g 4 + 15n 2 q 2 + 15) + O (M 4 )) , 

(4.38) 
with 

ql = q 2 a cS , a eS « a + 24C(3)a 2 M 2 . (4.39) 

We can read an 0(a 2 M 2 ) correction whose coefficient involves £(3) and C(5). 
For large M we use the following formulas 

Pl (M) = £(M) + 0(M- 2 ), £(M) = 21og(Me 7+1 ) , 

Pk (M) = 2((2k - 1) + 0(M- 2(fc - 1} ) , k > 1 . (4.40) 

Thus we find 

(e 2 ^ a ) = e" 2g * (1 + 27r 2 g 2 + O (M" 2 )) . (4.41) 

Now 

— ps - - 81og(Me 1+7 ) , (4.42) 

Oeff « 

which can be interpreted as the expectation value of the circular Wilson loop for Af = 2 
SYM written in terms of the renormalized gauge coupling. 

It should be noted that in the formulas ( I4.38p . (I4.4ip the additional a 3 corrections are 
suppressed, either by M 4 -in the small M limit- or by M~ 2 -in the large M limit. There- 
fore these formulas capture the full a dependence if we restrict to the leading behavior in 
the power series expansion in the mass (or inverse mass) parameter. 

4.2.2 't Hooft loop 

The expectation value of a supersymmetric 't Hooft loop operator carrying magnetic 
charge, with support on the equator of S 4 , is given in eq. (7.66) of [9]. We set the 9 
parameter of the gauge coupling r = ■%- + ^- to and restrict to the zero instanton 
sector. The magnetic charges are B = idiag(p/2, —p/2) and, for concreteness, we shall 
consider the case of p = even. Using the relations 

T(z)G(z) = G(l + z) , zY[z) =Y(l + z) , (4.43) 

we find the following form for expectation value of this ' t Hooft loop operator: 

{Tip)) = (Z-^ 2 *)- 1 J2 



x 



(E=l)\ (E±l\\ 

q=p,p-2,...,-p V 2 / - V 2 / - 

.i*A£ f q (a,M) h p (a,M)H 2 (2a) 



J aa 6 H(2a - M)H(2a + M) ' [ ' 
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where 

'sinh(vr(2a + M)) sinh(vr(2a - M)) 



h P (a,M) < . 

sinli (27raj 

cosh(47ra) — cosh(27rM) 
2smh 2 (2vra) 



(4.45) 



/o(a) = 4a and, for q even, |g| > 2, 



(4a 2 -M 2 )^ fi ((A; + 2M) 2 + 4a 2 )^ i - fe ((fc-iM) 2 + 4a 2 )i i - fe 
In particular, 

/ 2 (o,M) 

/ 4 (o,Af) 



4a 2 1 + 4a 2 



4a 2 - M 2 ' 

(4a 2 ) 2 (4a 2 + l) 2 (4a 2 + 4) 

(4a 2 + (1 - iM) 2 ) (4a 2 - M 2 ) 2 (4a 2 + (1 + iMf) 



These expressions show that for any finite magnetic charge p the structure of perturbation 
series is the same, due to the fact that poles appear only at imaginary values of a, 
not on the real axeso As in previous examples, the associated Borel series obtained 
by expanding the factor in the integrand multiplying exp[— a 2 /a] has finite radius of 
convergence, determined by the singularities of the integrand in the complex plane. The 
integration over a on each term of this series produces an additional factor n\. The 
perturbative series for (T(p)), which is again an alternate series, therefore diverges like n\ 
but it is Borel summable because the integral over a does not hit any singularity. 

For M = 0, the expression for the 't Hooft loop operator reduces to the one of M = 4 
SYM theory, studied in [22|[26H28] . Setting M = 0, we see that h p — > 1 and f q simplifies: 



2 



/ 9 (a,M = 0) = ^+4a 2 . (4.47) 



We get (p = 2m) 

«/■ - i * -V / 9m \ 4ir 2 n 2 f 167T 2 a 2 . „ 

(T( P ))m=o = {Z*=r i Yl {m + n) e ^j dae ~^( n +4a ) 



r2^2 






2 Note that the poles at 2a = ±M in f q (a, M) are canceled by zeros of h p (a, M). This special point is 
related to the point identified in [7], where a component of the adjoint hypermultiplet is massless. 
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Under an S duality transformation r — > — 1/r, 't Hooft loops are mapped to Wilson loops. 
The precise relationship is discussed in [D]. For 6 = 0, the S duality transformation is 
g — > 167c 2 /g and a given term in (14.48|) is mapped to the TV = 4 result (I3.27P with £ = n. 

The expectation value of the 't Hooft loop operator determines the behavior of the 
monopole-antimonopole potential. The large loop length behavior of the 't Hooft loop 
operator is believed to indicate whether the theory is in a confining, Higgs or Coulomb 
phase. For a confining phase, the 't Hooft loop operator is expected to exhibit a perimeter 
law, whereas for the Higgs phase it should exhibit an area law. In the present context the 
't Hooft loop is placed around the equator of S* 4 and it has therefore a length proportional 
to r, the S 4 radius. Although the scale of the 't Hooft loop cannot be disentangled from 
the scale of S 4 , it is nevertheless interesting to examine its behavior as r is changed with 
respect to the mass parameter. In this manner by large 't Hooft loop one means r 3> 1/M 
or, since we have set r = 1, it means looking at the behavior at M ^> 1. 

Let us first expand (T(p)) in powers of g (here we only study the g <C 1 regime, since 
we are ignoring instanton contributions). At g -C 1, in the sum over q, the terms with 

2 2/2 

q = ±p dominate because of the factor e 71 " q ' 9 . In order to study the behavior at M 3> 1 
it is convenient to express the one-loop determinant factor in terms of B^ =2 *(a 2 ,M) 
defined in (I4.34p . In particular, for p = 2, we need to compute 



4tt~ 



2e~9 r f , . , 9n 4a 2 (cosh(47ra) — cosh(27rM)) ^_ A . A _ 9 , , 9 , ,. 

(T(P = 2)) = w ^J d a(l + ^) lJ ( J a)(ia2 _ { M . 2) e «« B^ (a\ M) 

(4.49) 
a e ff has been defined in (I4.35p . We find 

(T(p = 2)) = 2e^ (T (1) + cosh(27rM) T (2) ) , (4.50) 



with 

T (i) 

T (2) 



1 3+(3 + 5tt 2 )M 2 



47r 2 M 2 a eff 6vr 2 M 4 

n " ,r (30M% + vr 2 (25 + vr 2 ) M 4 + 5 (3 + 5vr 2 ) M 2 + 15) + 0(a 2 ) 



5tt 2 M 6 

1 3-(tt 2 -3)M 



4vr 2 M 2 a cff 6vr 2 M 4 

+ ^^ (30M 4 p 2 + vr 2 (vr 2 - 5) M 4 - 5 (vr 2 - 3) M 2 + 15) + 0(a 2 ) . (4.51) 

The origin of the l/a e s behavior in the leading term can be traced back to the fact that 
the integrand in the numerator now contains a term 4a 2 + q 2 /A = 4a 2 + 1, whereas in the 
denominator Z^ =2 * contains the usual Vandermonde factor 4a 2 . 
In the large M limit we therefore find the following behavior: 

4^ 2 e 27vM a 2 

(T(p = 2)) = 4e^ — (1 -j- 2 log M + 0(g 4 ) + 0(M~ 2 )) . (4.52) 
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This is readily extended to general p. For the dominant q = ±p term at g <C 1 we now 
find 



i wprM / 2 \ 

(Tip)) = c(p) w —^ 1 - fa log(M) + 0(g A ) + 0(M- 2 ) , Mr » 1 

(- e -27rlm(r)j^4-)— <? V ^ / 

(4.53) 
where r = 4- + ^? — > ^B. and we have restored the r dependence. 

2tt g 2 g 1 r 

In [7], for the theory on R 4 , it was pointed out that the limit to pure N = 2 SYM 
theory is achieved from the N = 2* SYM theory by taking M large and holding e 2mT M 
fixed. Strikingly, in (J4.53P exp(27ur) appears in the right combination with the power 
of M. We now see that in our context the prescription by Seiberg and Witten has a 
simple physical explanation. The coupling constant dependence must be through the 
renormalized gauge coupling, 

ir?( 1- l? ,og(M) )- (4 - 54) 

At large M, there should not be any other dependence on the bare gauge coupling which 
is not through the combination forming g^. Indeed, (J4.53P can be written as 

i p -KprM 

{Tip)) = c{p) -, —^- (1 + O(g') + 0(M- 2 )) , Mr » 1 , (4.55) 

( e -27rIm(-r R n— 9r 

The remaining dependence on M is in the exponential factor e KprM . This originates from 
h p , which in turn comes from two contributions, the one- loop determinant at the equator 
and a contribution from monopole field configurations that screen the monopole charge 
B. Viewing N = 2* SYM theory as pure N = 2 SYM with a UV cutoff M, this factor 
might be indicating a very singular behavior of the expectation value of the 't Hooft loop 
operator in the UV limit, or it might be just reflecting the fact that the loop becomes 
very long and it follows a perimeter law. To clarify the meaning of this behavior it would 
be useful to compute the expectation value of a 't Hooft loop operator with support on a 
circle of arbitrary length L. This would allow one to consider a limit Mr ^> 1 at fixed L. 

5 Concluding remarks 

One of the lessons of the matrix model examples studied here is that the divergence of the 
perturbation series is not inherent to the theory but rather produced by an illegitimate step 
in the Feynman-diagram evaluation of functional integrals, when commuting an integral 
with a sum: the power expansion of vertex operators with the functional integration. In 
this respect, the Borel summation emerges not as an artificial trick "to make convergent 
what is actually divergent" , but it is rather built in the theory. In the simplest SU{2) cases 
the exact gauge theory partition function automatically appears in the form of a Borel 
transform. It is tempting to conjecture that gauge theories are free from divergent sum 
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pathologies if care is taken when commuting the functional integral with other integrals or 
infinite sumsjfl More precisely, for certain observables the Borel function may have poles on 
the real axes. In this case the Borel technique requires deforming the integration contour, 
a procedure that is not free from ambiguities and typically produces imaginary parts. 
Thus, a more precise conjecture is that there should be no divergent sum pathologies in 
observables where an imaginary part is not expected. For observables where an imaginary 
part is expected (like e.g. two-point functions of unstable particles) a suitable contour 
deformation prescription may not only render the integral convergent but also give rise 
to the expected imaginary part. 

To conclude, we summarize the main results. 

• The perturbation series for the partition function in the U(N) x U(N) ABJM the- 
ory is an alternate asymptotic series with (— l) m (2m)\(2/irk) 2m asymptotic behavior. 
This behavior is not associated with instantons (which are absent in three dimen- 
sions) but with poles on the complex plane of the one-loop determinant. The integral 
defining the partition function is convergent, so the divergent perturbation series is 
Borel summable. For the U(2) x U(2) theory, a simple integral representation illus- 
trates the origin of the divergent perturbation series. 

• For Af = 2 577(2) SYM theory with four massless hypermultiplets, we computed 
the asymptotic behavior of the perturbation series for the partition function and 
Wilson loop operator, in the zero and one instanton sector, and again found a Borel 
summable perturbation series, with asymptotic behavior (— l) n n\a n , a = g 2 /lQir 2 . 

• For the Af = 2* theory in the planar limit, the Wilson loop at weak coupling has 
the same form as in Af = 4 SYM theory upon replacing the 't Hooft coupling A by 

* /.^ x A \K'(M)\ XK'(M) 

MM) = A (l + -£->■) , -^ « 1 . 

At strong coupling and small mass parameter, we find W(C C i rc \ e ) ~ e 2?rT with 



V\ / M 2 



where T can be viewed as the tension of a semiclassical string configuration. In all 
cases augmenting the mass parameter produces an increase of W(C c i TC i e ). Finally, 
we have examined the expectation value (T(p)) of a circular 't Hooft loop operator 
of magnetic charge p/2 in the SU(2) theory. Like for the Wilson loop, the asso- 
ciated perturbation series in the zero instanton sector is divergent with a (— l) n n\ 
asymptotic behavior, and Borel summable. At large masses, where the theory is 
expected to flow to pure Af = 2 SU{2) gauge theory with a UV cutoff M, (T(p)) 
seems to increase as e nprM in the weak coupling regime. 



3 For example, in the U(2) x U(2) ABJM model, one can compute the integral (|2.26p by expanding 
S(v/k) in power series in v/k -which would lead to a divergent series - or by residues |18) - which would 
lead to a finite sum of terms. 
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It would be interesting to study sectors of higher instanton charge in M = 2 SU(2) 
SYM theories to see if more complicated instanton configurations may give rise to singu- 
larities in the Borel function lying on the integration region. This could provide an arena 
to clarify issues of ambiguities in contour deformation in the context of non-abelian gauge 
theories from closed analytic expressions. 

It would also be interesting to compute observables in N = 2* SU(2) SYM theory with 
a characteristic length L. By introducing a scale one can separate UV and IR physics in 
the series expansions and, in particular, see if the use of perturbation theory is problematic 
in the IR region, like in the case of QCD renormalons. 
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